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Iterated conformal mappings are used to obtain exact multifractal spectra of the harmonic measure 
for arbitrary Laplacian random walks in two dimensions. Separate spectra are found to describe 
scaling of the growth measure in time, of the measure near the growth tip, and of the measure away 
from the growth tip. The spectra away from the tip coincide with those of conformally invariant 
equilibrium systems with arbitrary central charge c < 1, with c related to the particular walk chosen, 
while the scaling in time and near the tip cannot be obtained from the equilibrium properties. 



Diffusion-limited aggregation jy and the dielectric 
breakdown model [|| (DBM) are central problems in frac- 
tal growth which still pose difficult problems. As the 
growth is controlled by a Laplacian field, the harmonic 
measure on these clusters is closely connected to the dy- 
namics. Quantitative characterization of the measure is 
provided by the multifractal exponents || . 

There are unfortunately only a few exact results for 
these exponents. The Makarov scaling law Q is true for 
any simply connected curve, and the electrostatic scaling 
law || is exact for the DBM. In an important recent 
advance, exact exponents were found for systems that 
may be described using conformal field theory |^](CFT); 
a series of spectra were found, labelled by the central 
charge c of the theory. This includes both equilibrium 
statistical mechanics clusters, such as Ising clusters at 
c = 1/2, as well as certain growth processes, such as 
random walks at c = 0. 

However, these exponents for growth processes only 
describe the static properties of these systems. Inter- 
esting dynamical questions, such as fluctuations about 
the average growth rate, remain unanswered in the ab- 
sence of exact results. To address this, we will consider 
a noncquilibrium growth process, the Laplacian random 
walk JjJ (LRW) , and show that its static properties away 
from the growth tip coincide with those of CFT clus- 
ters, but we will also obtain nontrivial scaling near the 
growth tip, giving rise to nontrivial dynamical scaling. 
This should shed light on more complicated growth pro- 
cesses such as DLA where the statics and dynamics are 
so deeply interwoven. 

The LRW is defined by a growing walk on a lattice, in 
which at every stage of growth a Laplacian field is com- 
puted with boundary conditions such that it vanishes on 
the walk and grows logarithmically at infinity. From this 
field an electric field E is computed, and the probability 
of selecting any site near the tip for growth is propor- 
tional to E v at that site, r] > 0. This can be considered 
as a dielectric breakdown model in which growth occurs 
only at the tip. 

In this paper, we propose a model of iterated confor- 
mal maps for the LRW. This model is not conformally 
invariant at the microscopic level, having a fixed parti- 



cle size; we then show that at large scales this model 
is equivalent to another discrete model which possesses 
microscopic conformal invariance. We obtain the exact 
multifractal spectrum as a continuously varying function 
of a parameter m, which will be indirectly connected to 
the r\ parameter above. The rj = I LRW corresponds to 
m = 1/2 or a c = -2 CFT [§. For other LRWs, there 
is a continuous relation between m and c for exponents 
of the measure away from the growth tip. While our 
model has the same continuum limit as the stochastic 
Loewner equation we use the discreteness to obtain 
additional exponents for scaling in time and measure near 
the growth tip. 

Conformal Mappings for the LRW — We use the 
method of iterated conformal mappings |l(| to construct 
an off-lattice LRW argued to be in the same universality 
class as that defined above. In the LRW, average growth 
follows the direction of the electric field at the tip of the 
walk, but with fluctuations about this preferred direc- 
tion. We will construct a model with the same features. 

Let F(z) be a map from the real line to a growing 
walk, with F analytic in the upper half plane. Define 
fx,\{z) to produce a bump of linear size y/~A at point 
x on the real line. For example, /o,a = V z 2 — A. If 
F(x) is the tip of the walk, then F(f x ^\(z)) maps to 
the walk grown by a distance \F'(z)\V\ in the direction 
of the field at the tip. Thus we propose the following 
model: fix the growth tip to be the image of z = 0. At 
each growth step, first compose F{z) — + F(f _\(z)), to 
grow the walk, where A = |F'(0)|~ 2 ^ with l\ a length 
that determines the size of the growth step in physical 
space. Then compose F(z) -> F(z±v r \l 2 /h), where the 
plus or minus sign is chosen randomly. This will shift 
the growth tip by a distance l 2 away from the preferred 
growth direction; the greater the ratio ^2/^1 the higher 
will be the dimension of the random walk and so the 
smaller the effective rj. 

We wish to determine the average of the q-th moment 
of at a point z = x + iy; that is, (\F'(z)\ q ), where 

the brackets denote ensemble averaging. This describes 
the field as a function of cutoff in the mathematical plane, 
z, which we we will later convert to scaling in the phys- 
ical plane. Denote this average after n growth steps as 
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Pq(x,y). By the chain rule we have 



We define r y (q) by 
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where we have used the approximation f(z) = z-X/(2z), Define a Legendre transform f y {a y ) = qa y - T y (q), with 



valid at scales larger than l\. 

Taking a continuum limit in the number of growth 
steps gives the differential equation 



dp q {x,y) _ A 



dn 



y 



+q 



(x 2 + y 2 ) 2 



(2) 



where we have replaced the discrete translation x — > x ± 
"vXfo/h by a diffusion term with coefficient m = (h/h) 2 , 
which defines the "meander" of the walk. 



Recall that A 
dt = Xdn in cq 



F'(0)\~ 2 lf. However, by rescaling 
), the dependence of A on |.F'(0)| 
can be removed so that we may assume A is constant. 
This rescaling does not alter the steady state solution 
of eq. (||). For growth in radial geometry, the absorber 
radius of the walk is e* = F\ , while in a cylindrical ge- 
ometry, the height of the walk is t. Once A is held fixed, 
eq. (^|) is the continuum limit of a growth process in which 
at each growth step F(z) — > F(fo tSl (z ± S2)), so that the 
fixed cutoff in physical space is replaced by a fixed cutoff 
si in the mathematical plane. This yields a construc- 
tion of the growth process as a random combination of 
conformal maps without memory (all maps chosen inde- 
pendently from one of two possibilities). This dynamics 
has a microscopic conformal invariance, demonstrating 
the conformal invariance of the LRW, which could have 
been anticipated by the known conformal invariance of 
the 77 = 1 LRW on-lattice. An image of a walk in cylin- 
drical geometry produced using this microscopically con- 
formally invariant dynamics is shown in Fig. 1, where the 
density of dots reveals the variation of the physical space 
cutoff. 

The ansatz p q = y^(l + x 2 /y 2 )^ solves eq. (§) 
with 



7i(9) =3 - 2m7 2 (g), 
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where we have chosen the sign of the square-root which 
gives the dominant contribution. 

Multifractal Spectrum — From the previous section, 
p q (0,y)(x,y) scales as near the growth tip (x = 

0), and as y 71 ^)- 2 72(«) awa y from the growth tip. To 
derive the multifractal spectrum r{q) in physical space 
from these, we define a multifractal spectrum r y (q) in 
the mathematical plane as a preliminary step, and then 
derive a useful inversion relation between T y (q) and r(q). 



T y (q). Then dividing the real line into segments 
of length y, the number of such segments which map 
under F into points separated by distance I — y ay in the 
physical plane scales as y~fy( a y\ Conversely, the number 
of boxes of size I which map, under the inverse of F, into 
size y — l a , scales as l~^ a \ Equating these yields the 
inversion relation for /(a) 



a=l/a y , f(a)=af y (a y ). 



(6) 



From eq. (^|) the functional inversion relation for r(q) 
follows, 



Ty(-T(q)). 



(7) 



We can break the integral in eq. (|5|) into a contribution 
away from the growth tip for x » y, and a contribution 
due to the growth tip for x < V- The first contribution 
alone would give T y (q) = 71 (q) — 272 (q) + q — 1, while the 
second alone would give r y (q) = r ji(q) + q- While the true 
T y (q) is obtained by the minimum of these two values, we 
choose to consider these contributions separately, giving 
moments away from the tip and near the tip. 

The Legendre transform of r y (q) away from the tip is 
f y (a y ) = 1 - (1 + m) 2 [(2 - a y ) + 1/(2 - a y ) - 2]/(2m). 
This is maximum at a y = 1 (Makarov's theorem pi), 
with / y (l) = 1 (Gauss's law). The inversion formula (|) 
gives 



/away (a) = OL + 
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This coincides with the CFT spectrum for c — 13 — 6m — 
6/to. Thus, c = 1 implies m = 1. For self-avoiding walks 
(c = 0) m = 2/3 and for 77 = 1 on-lattice LRWs (c = -2) 
77i = 1/2. For m < 1, the fractal dimension D = l + m/2, 
so the maximum dimension is 3/2. 

The Legendre transform of T y (q) near the tip is 
f y (a y ) = l + m-{2-a y )(l+m) 2 /(2m)-(m/2)[l/(2- 
cty)]. In this case, the inversion formula (p) gives 



/tip (a) 
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2 + 771 (2 + 77i) 2 (2a-l) , . 
+ — -- 8m (9) 



This has a maximum of zero at a = (1 + 771)/ (2 + 771). We 
find the superuniversal result that for a < 2/3, f(a) is 
dominated by the tip contribution (|^), while for a > 2/3, 
/(a) is dominated by contributions away from the tip. 

Finally, for fixed y << x, p q (x,y) oc a; 272 . This gives a 
further set of exponents for decay of the field along the 
cluster as a function of distance x from the tip at small 
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y, which differ from the exponents for decay of the field 
at x = as a function of y. In this case a Legendre 
transform of 272(g) yields a probability 

a ,-a(m+l) 2 /2m-l/(2ma) + (m+l)/m qqs 

of having a field of order £~ Q , with a > 0. Note that 
while the most likely value is a = l/(m+l), it is possible 
to have fields which are of the same order as the field at 
the tip. However, there is no scaling relation to connect 
these exponents to the decay of the field as a function 
of distance from the tip in physical space at fixed small 
distance away from the cluster. The reason we had the 
inversion relations (Q) above is that at a distance y from 
the cluster, the field is smooth on a scale y; however, for 
y << x, the field is not smooth on the scale of x. 

The exponents (^) are invariant under m — > 1/m, al- 
though the other sets of exponents are not. For m > 1 
the diffusion term is sufficiently strong that the growth 
process does not give rise to linear walks, but rather clus- 
ters, as a point with j/«0 and large x can diffuse back to 
x = despite the advection term in eq. (^) (see Figs. 2,3). 
While the dimension of such clusters increases above 3/2 
as m increases above 1, the dimension of the perimeter 
decreases as m increases and remains below 3/2, giving 
the invariance of exponents away from the tip mentioned 
above. Unsurprisingly, then, that the exponents near the 
tip differ and depend on m and not just on c. 

Dynamical Scaling — There is a form of the elec- 
trostatic scaling law || for this system. We have 
(in radial geometry) {dn/dF^ = l/(A J F 1 i)«. Tak- 
ing averages on both sides, and noting that (A~ 9 ) oc 
{F l /l 1 ) 2 i{h/F 1 y(~ 2 ^+ 2 <i oc F- T[ - 2q \ we get 

{{dn/dFxY) ocF- q ~ T{ ~ 2q \ (11) 

where r(q) is the scaling near the tip, obtained from Leg- 
endre transform of eq. (||). Now consider the average 
mass, n, at given radius, F\, Integrating eq. (fill), we ob- 
tain (n) oc F ± T( - . Since — r(— 2) = l + m/2, form < 2, 
the dimension obtained from the electrostatic scaling law 
agrees with that obtained above. This also gives a dimen- 
sion greater than 3/2 for clusters with m > 1. However, 
eq. ( |ll|) shows that there are large fluctuations in the 
growth rate, and that the typical growth is slower than 
average. 

Now we consider the time-averaged growth rate. Af- 
ter transforming to time t, eq. (Q) has only simple scal- 
ing in the mathematical plane, t oc ^fy. As a first 
step to finding moments of the time-averaged growth 
rate, consider the average of a product of \F'\ at dif- 
ferent times: (\Y\yiF'(x — 0, yt, ti))\ q ). Suppose all 

i 

|ij — tj\ are of order y 2 . » yf, so that the derivatives 
are correlated only down to scale yo, and uncorrelated 
below. To determine the effect of the correlations, con- 
sider this average for a single time but suppose F'(0,yo) 



is given, which imposes boundary conditions on eq. (|2). 
Then, {\( V F\Q,y)\9) oc (y/yo) T * (q} \F'(0, y )\ q . Thus, for 
r points at separate times the ensemble average is 

(HI^'^'W' *••))!*) ^Y[(yM Tyiq) (\ F '^yo)\ qr )- (12) 

i i 

Transforming to the physical plane we find (given that 
there is probability {yi/yo)~^ v ^ av ^ of y t mapping to U oc 
(yi/yo) ay \F'{0, yo)\ and using arguments as above eq. (||) 

i i 

We can now obtain moments of the change in mass 
(ri2 — n\) q between times ti,ta by integrating J dn/dt 
and applying eq. (|l3|). 

(( Jdn/dty) oc F- qT{ - 2 \t 2 t^-*^ 1 . (14) 

ti 

t 2 

For h-ti> 1, this is ((/ dn/dty) oc F~ qT( ~ 2 \ so that 
*i 

the total mass at given radius is self-averaging. Over 
shorter times, eq. ( [l4] ) indicates multiscaling of the mass- 
radius relation, or dynamical multiscaling connecting the 
two times t, n. The origin of these large fluctuations is 
that the tip may turn and grow towards the center of the 
cluster rather than away, so that F\ may increase very 
slowly with mass. The fluctuations in F\ are much less 
in DLA jll]], where growth occurs over the whole surface; 
while \F' | -2 may fluctuate locally, the average of |-F'| -2 
over the cluster has much lower fluctuations. 

Conclusion and Renormalization Group — We have 
demonstrated the conformal invariance of the LRW and 
calculated exactly the multifractal exponents. This pro- 
vides an opportunity to realize clusters of different con- 
formal field theories as different LRWs, giving an inter- 
esting connection between the iterated conformal map 
technique and CFT. We obtained scaling of harmonic 
measure both along the curve and near the growth tip. 
Because we have a solvable dynamic model, we are able 
to obtain exponents characterizing the scaling in time. 

After rescaling of time, the growth process was de- 
scribed by random composition of two different confor- 
mal maps. Other fractal sets such as the Julia set are 
obtained by iterating a single conformal map. 

It is interesting to note the lack of universality in these 
results, so that by changing the ratio h/h we are able to 
obtain continuously varying exponents. In the on-latticc 
LRW defined above, if successive particles are allowed, 
instead of attaching to the tip, to attach within some 
distance from the tip, we expect the dimension to change 
as the distance is varied, even with 77 held fixed. 

Recently, it was arg ued Q that for 77 > 4 DBM 
clusters are asymptotically branchless, and hence one- 
dimensional. One may wonder if instead these branchless 
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clusters are LRWs, as suggested previously Q, and have 
dimension greater than one. However, the lack of univer- 
sality in the LRW implies that even if such DBM clusters 
with given 77 were to behave as LRWs, the resulting di- 
mension need not be the same as that of an on-lattice 
LRW defined with the given 77; this 77 may be renormal- 
ized. More importantly, we claim that an asymptotically 
branchless DBM cannot behave as an LRW; if such a 
cluster were to behave as an LRW, then the tip would 
occasionally double back on itself, so that the point of 
strongest electric field would not be at the tip (as shown, 
the field at large x in an LRW can be of order the field 
at the tip), causing a new side branch to appear and 
contradicting the assumption of no branching. Thus, the 
possibility of side branching renormalizes the meandering 
of the LRW to zero and what results is a single straight 
branch. The asymptotic approach of the dimension to 
one for DBM clusters with 77 > 4 is confirmed numeri- 
cally @. 

Formally, one may extend the RG by noting that the 
change in growth trajectory induced by a single tip- 
splitting event gives a logarithmic correction to the clus- 
ter size, an additional contribution beyond the correction 
proportional to branch size considered before p2| . This 
yields, at lowest order in tip-splitting, a log-squared cor- 
rection to the cluster size, which is accounted for by an 
additional scaling field, representing the effective m of 
the LRW described by the cluster. However, as argued, 
the possibility of side branching renormalizes the m to 
zero. Thus, we may ignore these log-squared corrections, 
as done previously. 
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FIG. 1. Plot of an LRW with m < 1, using conformally in- 
variant microscopic dynamics and ~ 50, 000 particles. Points 
indicate position of tip after each growth step. 




FIG. 2. Plot of an LRW with m > 1, using fixed cutoff in 
physical space and ~ 50, 000 particles. 
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FIG. 3. Plot of an LRW with m > 1, using conformally 
invariant microscopic dynamics and ~ 50, 000 particles, m is 
greater than for the cluster in Fig. 2. 
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